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We discuss the effects of magnetic impurities on d-wave superconductors. We calculate the electron 
mediated RKKY interaction between the impurity spins in a d-wave superconductor and find that it 
decays as at large distances. We argue that this interaction leads to the formation of a spin glass 
at low temperature T <^ Tc- It was previously shown that a local complex ~ dxy order parameter 
is induced around each impurity spin. We consider the pair tunneling resulting in the Josephson 
interaction between different patches of induced dxy order parameter. Due to the local coupling 
between impurity spins and the superconducting order parameter the Josephson coupling favors a 
ferromagnetic phase at low temperatures. The competition between the Josephson coupling and the 
RKKY interaction gives rise to an interesting phase diagram. At low impurity concentrations we 
find an unusual supercondcting phase glass, where the impurity spins and dxy component are 
disordered and yet the product of these two develops a true long range order (S'^A^). This phase 
has no analog in purely magnetic spin glasses and arises as a result of the direct coupling of the 
impurity spin to the phase of A^. At high impurity concentrations it is possible that a ferromagnetic 
phase will form. 

PACS numbers: 74.25.Bt, 74.62.Dh, 75.10.Jm 



I. INTRODUCTION 

It has been known for a long time that magnetic ion^ 
in a superconductor can form a spin glass (SG) state.EJ 
There is growing experimental evidence for the formation 
of a SG phase in the high-Tc materials once magnetic, 
ions, such as Fe and Ni, are doped into the system.El 
Since magnetic ions interact witli. the superconducting 
condensate it is natural to expectu that this interaction 
will cause frustration of the underlying condensate and 
eventualy might lead to a superconducting glass or phase 
glass (PG). 

In all of the above discussions of the role of impurity 
spins it was assumed that magnetic scattering frustrates 
and suppresses the superconductivity. There are addi- 
tional physical effects that were not addressed in previ- 
ous work: namely how frustrated localized spins in an un- 
conventional superconductors can distort the condenstate 
and produce patches of secondary components of the or- 
der parameter near the localized impurity, see Fig. |l|. It 
has been shown that a magnetic impurity in a dx^-y^ 
superconductor with order parameter A° induces a local 
complex dxy conmonent of the order parameter which 
we designate A^.a Locally near each impurity site, on 
the scale of the coherence length ^, there is a patch of 
A'^ -I- iA^ order paramter for = -1-1 and A'^ — iA^ 
for = —1, depending on the sign of the impurity 
spin. Hereafter we assume classical spins with S — 1. 
This is a reasonable assumption taken the fact that mag- 
netic ions substituted into high-Tc superconductors have 
a large spins, such as Ni (S=l). 

Here we will explore the coupling between different 
patches due to the Josephson coupling. Phase coherence 
between these patches of A'' ± lA^ would lower the ki- 
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FIG. 1. Three impurity spins at sites i, j and fc in a 
two-dimensional dx2_y2 superconductor. The patches around 
site j and k indicate the induced d' — dxy component of the 
order parameter. Note that there is no induced order param- 
eter around spin i, since the spin is pointing in the xy plane. 



netic energy of the condensate and would tend to align 
all the patches into globally coherent A°-|-iA^, or its con- 
jugate. This would imply ferromagnetic ordering of the 
impurity spins . On the other hand, the dipolar and 
RKKY spin-spin interaction terms, which we show are 
mainly antiferromagnetic, would frustrate this ferromag- 
netic order. In fact, we find that the spin-spin interac- 
tion is frustrating and produces a SG phase at T < Tsg 
in the absence of a coupling between the superconduct- 
ing condenstate and spin degrees of freedom. When this 
coupling is considered the result depends on the relative 
strength of the two interactions. We find that since the 
spin-dependent part of the Josephson coupling is small 
the magnetic susbsystem will drive the phase transition, 
at least at low and intermediate impurity concentrations. 
Given that the impurity spins form a SG phase we are led 
to the question how the phase of the induced component 
±zA^ is affected by spin frustration. We find that the 
spin frustration in the SG phase will frustrate the phase 
of the dxy component and a PG will form at low tempera- 
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tures T < TpQ < Tsg- More generally we will discuss the 
possible phases of the coupled spin-phase model, and at 
higher impurity concentrations we will consider the for- 
mation of a ferromagnetic phase. When discussing glassy 
phases, we will consider states where the spatial average 
(5^) = 0. However at any given patch, i, the time aver- 
aged (S"? )r ^ 0. Similarly we will consider site and time 
averaged induced component (A^) = and (A^^),- ^ 
respectively. 

The additional new order parameter we find relevant 
is the product (S'^A^). The new order in this case arises 
from the possibilty of having fully disordered phases 
(A^) = (5'^) = and still having the true long range 
order in {S^A^) ^ 0. Physically it corresponds to phase 
locking of the patches to the local value of impurity spin 
S^. This order parameter has no analog in purely mag- 
netic spin glass systems and is a direct consequence of the 
spin-phase order parameter coupling discussed above. In 
all of the discussion hereafter we assume that the "back- 
bone" order parameter A° has true long range order as 
it is robust at low temperatures T <^Tc ^ lOOK. 

The plan of this paper is as follows: in the first section 
we explore the effective spin-spin coupling and integrate 
out the superconducting degrees of freedom by perform- 
ing an RKYY calculation. The effective spin model is 
found to give rise to a SG phase. In the second section 
of the paper the Josephson coupling between different 
patches is considered. It is found that the interaction 
favors a FM phase. In the last section we dicuss the 
PG phase, where the induced superconducting order pa- 
rameter is locked to the impurity spin. Furthermore we 
will discuss the possibilities of para- and ferromagnetic 
phases. _ 

In related recent workH Simon and Varma treat the 
magnetic impurity problem by a variational approach. 
They concentrate on the single impurity, but conclude 
by arguing that the formation of a ferromagnetic phase 
is unlikely and that a spin glass state is much more likely. 

II. EFFECTIVE SPIN MODEL 

In this section we will consider the effective spin model 
that describes the interaction between impurity spins sur- 
rounded by patches of induced complex dxy order param- 
eter, as shown in Fig. The most relevant sources of 
interaction between the impurity spins are the electron- 
mediated RKKY-interaction and the direct dipolar niag- 
netic interaction. The standard RKKY interactionO is 
mediated through an interaction of the form JS ■ /, where 
S referes to the spin of the conduction electron, and / 
denotes the spin of the impurity. Of particular interest 
to this work, however, is another term that is directly 
related to the induced order parameter. The L •/ in- 
teraction, where L refers to the angular momentum of 
a conduction electron, scatters an electron into the com- 
plex dxy phase. A second order RKKY calculation where 



the interaction potential is taken to be of the from L ■ I 
will therefore be relevant in this context. The effective 
spin Hamiltonian is thus of the form 

H = H''^ + H^-' + H^', (1) 

where the dipolar term H^^ and the electron mediated 
interactions H^'^ and iJ^'^ are given by 

= E it^' • - 3(5^ • f^j)iSj ■ hj)] (2) 
^"■'-E4'^^-^. (3) 
H'^' -Y.4''StS^^ (4) 

The size of the dipolar term for two spins separated 
by 1 Acan be estimated to be about 10~'* eV, while the 
RKKY interaction depends strongly on the coupling be- 
tween the conduction electrons and the impurity spin. 
For Mn ions in alloys is has been estimated to be on the 
order of 10^^ eVJj We are not aware of direct estimates 
of the RKKY interaction in high-Tc materials, but it is 
likely to be greater than the dipolar interaction. The spa- 
tial decay of the last two terms are explicitly calculated 
in the appendix and are of the form 

^^.j ^ 2J^Iiqkl ^ ^ sin(2if r)] (5) 
^s-i _ J'^h -h jQ_3^ _|_ Q 33 sin(2Xr)] (6) 

where v is the gap velocity and K is a momentum cut- 
off. The RKKY terms thus share a cubic decay with the 
dipolar interaction, but are found to be completely an- 
tiferromagnetic. However, independently of the specific 
form of the terms, general consideration tells us that as 
long as there are sufficient amounts of disorder_£ind frus- 
tration present the SG phase should be realizedB The im- 
purity spins we are considering are randomly distributed, 
and antiferromagnetic interactions will therefore lead to 
frustration. So as long as the effective interactions are 
not primarily ferromagnetic the model should have a SG 
groundstate. The dipolar interaction is antiferromagnetic 
in the z-component of the spins in an xy plane, and ferro- 
magnetic in the plane. As mentioned above the RKKY 
terms are completely antiferromagnetic. We therefore 
conclude that the effective spin model, consisting of the 
combined dipolar and RKKY terms, exhibits a spin glass 
phase at low temperatures. For large enough quantum 
fluctuations it would also be possible to realize a quan- 
tum paramagnet phase, where any freezing is destroyed 
by quantum fluctuations. We focus here on large classical 
spins and assume that quantum fluctuations are negligi- 
ble. In principle similar considerations can be given for 
S" = 1/2 impurity spins, in which case the paramagnetic 
phase could be realized. 
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Note in particular that the impurity spins would form 
a spin glass even without the H'"'^ term, and that this 
term is second order in <C A° . This is the motiva- 
tion why we have first considered the effective spin model 
and determined the spin configuration. So although the 
spin and superconducting degrees of freedom are con- 
nected, as we will see in the next section, we will there 
take the spin configuration as given, and determine the 
phase of the superconducting order parameter based on 
this. Note also that the importance of the dipolar term 
may increase as the Ai component is induced, since this 
will open up a fully gapped state which will tend to ex- 
ponentially decrease the long-range RKKY interaction. 
The short-range RKKY interaction should not be signif- 
icantly affected by the opening of the superconducting 
gap. 

III. JOSEPHSON COUPLING 

In this section we will consider the Josephson coupling 
between the different patches of induced order parame- 
ter. Let us first consider the order parameter around an 
impurity spin 

= (AO + e'^^?Ai)e^^', (7) 

where the real parameter A'^ denotes the d^^-y^ order 
parameter, and the likewise real quantity A^ denotes the 
induced d^y component. Hereafter we assume that A° is 
a robust order parameter that develops true long-range 
order at T < Tc ~ 90 K and remains ordered in all of the 
phases we discuss. The impurity spin uniquely deter- 
mines the relative phase of the induced order parameter 
to be +TT/2 for S"^ = +1 and -tt/2 for 5^ = -1. The 
main result of the previous section was that the impurity 
spins, considered independently of the coupling to the in- 
duced order parameter, will exhibit a spin glass phase at 
low temperatures. Let us next consider what effects the 
inter-patch Josephson interaction may have. 

In order to address this question we examine the 
Josephson coupling between different patches given by 

H' = -j2\hjm^j + 'i'r^*), (8) 

Using the local order parameter Eq. (Q) this leads to an 
interaction of the form 

H' = -2j2mi^y + i^'yS^S]] cos{e., - e,) (9) 

We have here assumed that A^ = A^ = A°, and like- 
wise for Ai. The first term, which is zeroth order in Ai, 
wants to align the phase of A° at different patches by 
setting 9i = 9j . The second term wants to align the im- 
purity spin in an ferromagnetic phase. The Josephson 
coupling thus favors a ferromagnetic spin configuration, 
while the effective RKKY and dipolar spin model favors 



a spin glass. In the next section we will consider some 
possible outcomes of this competition. 

We would, however, like to point out that the RKKY 
and Josephson effects are not as independent as they may 
at first seem. Both are mediated by electrons and while 
the dominant part of the RKKY interaction in the super- 
conducting state is zeroth order in A^ there is a part that 
is second order in A^, corresponding to the Josephson 
coupling. The Josephson interaction physically expresses 
an effect arising from electron pair tunneling, while the 
GG part of the RKKY interaction expresses an electron- 
hole channel. Here G and F are the normal and anoma- 
lous propagators in superconducting state. The FF part 
of the RKKY interaction is closer to the Josephson in- 
teraction since it expresses an electron-electron process. 
Furthermore, there is a second order contribution to the 
RKKY interaction which we have not evaluated in this 
work. This contains two explicit d^y propagators, with 
two exchange interactions at the impurity spins. This 
part of the interaction should be similar to the Josephson 
interaction since explicit information about the patches, 
such as spatial decay, should be contained in the propa- 
gators. The physical effects of this term should, however, 
be included in the Josephson coupling considered above. 

IV. POSSIBLE PHASES 

In this section we will consider the different phases that 
could occur as a result of the competition between the 
RKKY and the Josephson terms. First we will enumerate 
the different phases and thereafter we will discuss them 
in some more detail. 

The simplest phase is a ferromagnetic (FM) phase, fa- 
vored by the Josephson coupling of different patches. The 
FM phase is characterized by a finite spatial average of 
the magnetization and the induced phase; (S^) ^ and 
(Ai)^O. 

The phase favored by the effective spin Hamiltonian, 
on the other hand, is a SG phase, where the spatial av- 
erage of the magnetization vanishes 

(5-) 0, (10) 

but the time averaged local magnetization remains finite; 

^Tlr + 0. (11) 

Assuming that the phase of A^ is determined by the im- 
purity spin this gives rise to similar ordering for the in- 
duced phase; (A^) = and (A^)^ ^ 0. In this case the 
combination S'^A^ will also have a non- vanishing spatial 
average 

(S'^Ai)^O (12) 

This particular order parameter is unique for the coupling 
of SG and SC degrees of freedom and is not present in a 
purely magnetic system. It describes the phase-locking of 
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the superconducting phase and the impurity spin.EI There 
also exists the possibility that the induced order param- 
eter does not phase-lock with the impurity spin, but 
prefers to varypO,ver length scales greater than the inter- 
patch distance.ll2l In that case the spatial average (S^A^) 
would vanish, but the local time average (S'f A^),- would 
still remain finite. 

If quantum fluctuations are strong enough, then it 
would also be possible to realize a paramagnetic spin 
phase. In this case also the time averages {S'')t and 
(A^)t- would vanish, and only local (S^A^) would remain 
finite, for as long as the phase-locking is maintained. In 
this work we will, however, neglect this possibility since 
we focus on large classical spins. 

After having considered the different options let us con- 
sider the interplay between the different terms. If there 
is a low concentration of impurity spins (the patches do 
not overlap), then the ferromagnetic Josephson term is 
bound to be exponentially small, and the ground state 
should be a spin glass, characterized by a no n- vanishing 
spatial average (S^A^). If the impurity concentration 
becomes larger, but not large enough to kill the A" su- 
perconductivity, then there may be a region where the 
ferromagnetic Josephson terms are predominant and the 
FM phase is formed. The FM phase was considejred in 
detail in a previous Ginzburg-Landau descriptionO. 
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FIG. 2. Phase diagram for the impurity doped d-wave su- 
perconductor as function of impurity doping (x) and temper- 
ature (T). The labeled diagrams show the normal (N) phase, 
the d^2-y2 superconducting (SC) phase, the spin-glass (SG) 
phase, the superconducting d^2_y2 + id^y phase-glass (PG) 
phase and the ferromagnetic (FM) phase. The SC phase is 
suppressed by disorder due to pair-breaking, while the SG 
phase is independent of the SC order parameter and will per- 
sist in the N phase. The PG phase is induced by the SC and 
SG phases and exists only within these phases. The possible 
FM phase is a result of strong Josephson coupling when the 
impurity concentration is large. It is strongly suppressed close 
to Xc due to the suppression of the superconducting order. 

In Fig. |2| we present a phase diagram as a function 
of impurity concentration and temperature. Note that 



this PG phase is different from previously proposed su- 
perconducting phase glasses in that the glassy behavior is 
only displayed in the induced component of the order pa- 
rameter, and not in the robust dx2^y2 part. Considering 
the phase diagram we note that the disorder suppresses 
the critical temperature for the superconductor-metal in- 
sulator, as is well known. The impurity spins form a 
spin-glass phase at low temperatures, and this phase is 
independent of the electronic order parameter and per- 
sists also in the metallic phase. The PG phase is induced 
by the SG and SC phases and hence it must only exist 
within the boundaries of these two phases. The possible 
FM phase is induced by a large impurity concentration. 
As the dj.2_y2 superconducting order parameter gets sup- 
pressed by disorder the induced component will also van- 
ish, and hence the FM phase gets strongly suppressed as 
we approach x^. 

Experimentally the proposed phase locked state can 
be observed in scanning tunneling microscope measure- 
ments, where the particle spectrum would develop a full 
gap near impurity site even though the phase of A^ re- 
mains uncertainQ. The ac magnetic susceptibility in the 
superconducting state also should show features upon 
crossing the SG and PG lines. Another experiment, 
that would be sensitive to the appearance of the (5'^Ai), 
would be the penetration depth that would become expo- 
nential 51/ \^ (X exp[-|Ai|/r] in the PG and FM phases. 

V. CONCLUSION 

We have examined the role of magnetic impurities in 
a d-wave superconductors. In particular we have stud- 
ied the effective impurity spin model arising from elec- 
tron mediated RKKY and magnetic dipolar terms and 
argued that these terms lead to a SG phase at low tem- 
peratures.. Furthermore we have analyzed the coupling 
between the spin- and superconductor order parameter 
that arises from the Josephson interaction of patches of 
induced order parameter around the impurity spins. The 
Josephson interaction favors a FM phase. At high impu- 
rity concentrations the FM phase may be realized, while 
at low concentrations the SG phase would be prefered. 
Due to the coupling between the spin and superconduct- 
ing order parameters the SG phase induces a supercon- 
ducting PG at low temperatures. The glassy behavior is 
a property of the induced d^y component of the order pa- 
rameter, while the primary component dx2-y2 is assumed 
robust. The superconducting PG phase is characterized 
by an order parameter of the form (S^A^), which de- 
scribes the phase-locking of the induced order parameter 
and the local impurity spin. In addition to the PG and 
FM phases we have also discussed a possible paramag- 
netic phases. 

We are grateful to D. Agterberg, N. Bonesteel, M. Graf 
and I. Martin, for useful discussions. 
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APPENDIX A: THE RKKY INTERACTION 

The RKKY interaction describes a second-order pro- 
cess, where an electron of momentum k interacts with an 
impurity spin at r = Ri . is scattered to a state with mo- 
mentum q and interac;ts witli a another impurity spin at 
r = i?2, where it is scattered back into the original state. 
This process is described by the following diagram: 
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FIG. 3. Second-order RKKY interaction. 

Using Feynman rules we get the following expression 
for the effective Hamiltonian 

d''q 



iH ■ 



d'^k r d'^q f duj 



(A2) 



{2ttY J {27tY J (27r) 



Using the Nambu formalism for the superconductor 
this expression transforms to 
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(A5) 
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where Ek = ye^-i-A|^. For a tight binding model Cfc = 
—t [cos{kxa) + cos(fcj^a)], and the rfx^-j/^ gap is given by 
Afc = Aq [cos{kxa) — cos{kya)]. 



Next we will consider a fcw specific forms of the in- 
teraction , which is the Fourier transform of the 
electron-impurity spin interaction. For free electrons 

\k) = e^' 



Vkl,,0 = = J d'^xe-'('=-«)-(a|F«|/3) (A9) 

Let us first consider a contact potential of the form = 
J6{x — R)S ■ I, where S denotes the electron spin and / 
the impurity spin. This results in 



14l,/3 = ^e-'('=-'')«(a|5.n 



(AlO) 



Next we will consider a, L-I interaction, where L is the 
angular momentum of the electron, and / is the spin of 
the impurity at location R. The electron moves relative 
to the impurity spin, which sees a magnetic field of the 



form B 



\x-R\ 



7, where x is the location of the electron. 



We will consider two spatial dimensions, and since the 
angular momentum of the electron will have only a z- 
component we consider an interaction of the form 



= J- 



L-I 

x-R\^ 



J- 
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(All) 



where L — {x — R) x p = —i{x — R) x V^. We get the 
following expression for the matrix element: 



7' 



re''"'{a\l3) 
(A12) 



After performing the integrals we arrive at 
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(A13) 



We are now in a position to evaluate the effective in- 
teraction 
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(2^0^ y (2^7 (2^ 
[Gii(w, k)Gii{LO, q) + G22(w, fc)G22(w, q) 
+2F{Lj,k)F{uj,q)]. (A14) 



We will begin with the frequency integral for the FF con- 
tribution 



If = 



duj 



2F{uj,k)F{oj,q). 



(A15) 



There arc poles at uj = i^yE'^ — iS = ^Ek ± i5- Closing 
the integral in the upper half plane leads to contributing 
poles at w = —Ek + i6 and co = —Eg + i6. It follows that 
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Next we will consider the GG contribution 



(A16) 
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Ig = 



duj 



Gil (w, fc)Gii (w, g) + G22 (w, fc)G22 (t^, <?) 



Proceeding as above it follows that 
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We will start by considering an interaction of the L cot / 
kind. We then have 
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where r = i?i — i?2- For the effective interaction we then 
get 
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In order to solve the abawe integral we introduce the 
nodal point approximation.EjThe dominant contribution 
to the integral should come from each nodal point, and 
we perform a rotation and translation to transform the 
origin to the nodal points, with the x-axis along the tight 
binding Fermi surface, see Fig. ^. 
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FIG. 4. Node notation. 



The transformation is given by 
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where kp = 7r/V2a and 6 = {7r/4, 37r/4, 57r/4, 77r/4} for 
the four nodes. 

Next we need to apply these transformations to all 
quantities in the effective interaction. The energy of a 
tight-binding model, tk = — t [cos(A:2;a) + cos(A:j,a)], will 
transform according to ej. — vpk'y for all nodes, and the 
gap function, = Aq [cos(fc;ra) — cos(fcyo)] will trans- 
form according to A^ — ±.V/\k'^, with a positive sign for 



nodes 1 and 3, and a negative sign for nodes 2 and 4. 
Assuming an isotropic superconductor v = vp = v/\ this 
leads to Ek = vk' . 

The rotation and translation cannot affect j/c — g|, and 
therefore |fc — = \k' — q'\. Furthermore 
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for all the four nodes. The second term is linear in the 
momentum and will be retained. The final term to be 
transformed, {k — q)r, is considered next. This term de- 
pends on the direction of r, and different nodes will give 
different contributions. Assume that r is fixed in a direc- 
tion 6r, and the contribution to the effective interaction 
from node 1 has been calculated. Since all other terms 
are identical, the contribution from the other nodes must 
be given by substituting {9r + tt/2, Or + 27r/2, Or -I- 37r/2} 
for Or in the expression obtained for node 1. Therefore 
we will look at how the transformation is done for node 
1, and the results for the other nodes will follow. For 
node 1 we get 



[k - q)r = rx{kx - fe) + ry[ky - qy) 

[cos Ok'-q' (cos Or + sill Or) 
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(A23) 



As a summary we have thus arrived at the following 
transformations 

efe = vk'y 
Afc = vk'^ 
Ek = vk' 
\k-q\ = \k'-q'\ 
{k X q)' = kp(k' -q')"" 
r\k' - q'\ 

{k — q)r = -j= — [cos 6*^/-^/ (cos + sin6'r) 

v2 

+ sinOk'-qi (— cos Or + sin6'r)] 

Dropping the primes, and using these results we can thus 
linearize the effective interaction around the nodes: 
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We start by integrating over the angles. First we fix the 
relative angle between k and q and integrate over Ok, 
thereafter we integrate over 0. Integrating over Ok we get 



Att^v J J k + qj 
TT[Jo{r\k - q\) - sm{20r)J2{r\k - q\)] 



(A25) 



The angular dependence will, however, vanish, since sum- 
ming up the contributions from the four nodes gives us 

sin(26'r) + sin(2(6'.r + tt/2)) 
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+ sin(2(6'^ + tt)) + sin(2(6lr + 37r/2))] = 0, (A26) 



and this tells us that the effective interaction is isotropic, 
even though the gap is anisotropic. Integrating out the 
relative angle we find 
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This integral is oscillatory, and we introduce a momen- 
tum cut-off K and make the integration variables dimen- 
sionless by letting k ^ rk and q ^ rq: 
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This integral can be solved numerically, and the behavior 
for large r is given by 
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This results represents an anti-ferromagnetic part 
that is independent of the cut-off, and an oscillating 
sin(2_fi'r)r~^ part. Due to the relative sizes of the two 
terms the interaction is always positive and licniec; com- 
pletely anti-ferromagnetic. The S ■ I interaction will lead 
to a very similar result, differing only in the prefactor. 
This can be seen, since the result, after performing the 
first angular integral will be given by 



H 



dk I dq- 



kq 



Att^v J J ^ k + q 
d0{cos9-l)'KJoir\k-q\ 



(A31) 



and the only difference compared to the L ■ I interaction 
is the prefactor. The final result for the S ■ I coupling 
will be 



H 



SI 



[0.37 0.33 sin(2ii:r)] 



(A32) 



So the two electron-spin couplings lead to the same func- 
tional form of the RKKY interaction. We have therefore 
showed that both the 5 • I and L ■ I interactions give rise 
to an effective anti-ferromagnetic model. We have used 
the nodal approximation and assumed that the supercon- 
ductor is isotropic vp = va, and these approximations 
may change the final result somewhat, but it appears un- 
likely that they would make the effective spin-spin inter- 
actions predominantly ferromagnetic, and therefore the 
randomly distributed spins will form a spin-glass phase, 
as discussed in the main part of the paper. 
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